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Fourier Coefficients of Mathieu Functions in 
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A method for calculating the Fourier coefficients of Mathieu functions in stable regions based 
on the method of bisection and Miller's recurrence algorithm is introduced. Some Fourier coefficients 
are calculated and compared with those given earlier by Tamir and Wang. It is shown that the method 
of Tamir and Wang fails for indices in the neighborhood of integers. 
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It is known from the theory of Mathieu functions [1] ' that the Floquet solution of the differ- 
ential equation 

y"(z) + (A-2/r cos 2z)y{z)=0 (1) 

with characteristic exponent v has the Fourier expansion 

oc 

me v {z,h z )= ^ c£ n (A 2 )e I " (l ' +2 ">*. (2) 

rc=-oc 

The Fourier coefficients satisfy the following three-term recurrence formula: 

-h*c v 2n+2 (h*) + [X- (i/+ 2/i) 2 ]cJ n (A 2 ) -A 2 e* n _ 2 (A 2 ) = 0. (3 ) 

Continued fraction expansions for the eigenvalues k=k(v, fi 1 ) and for the Fourier coefficients 
c 2 W / c o are g iven by Tamir and Wang in [2j. These relations are valid both in the stable (v real) and 
unstable (v complex) regions. F(v) = c%J<% is numerically evaluated and depicted. It is seen from 
the graphs in [2] that the given formulas for c»Jc% are not particularly suitable for numerical evalu- 
ation in the stable regions. Representations for F(v) = c\Jc\ derived from other numerical methods 
will now be introduced. 

The eigenvalues k of Mathieu's differential equation can be determined as the eigenvalues of 
the symmetrical tridiagonal matrix of the homogeneous system of eq (3) with the help of the bisection 
method [3|. Then the Fourier coefficients are obtainable from the recurrence formula (3) under 
consideration of stability (Miller's recurrence algorithm [4,5]). The accuracy of this method depends 
principally upon the number of equations used. 

Numerical work was carried out to single precision on the computer of the Technische Hoch- 
schule Darmstadt. For evaluating the eigenvalues X, a modified version of the bisection method [6] 
was used. This modification prohibits overflow. Also, a numerical test ensures that the trunca- 
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tion error is smaller than the rounding error of the bisection method, when the number of equa- 
tions used is [v + h 2 ]-\-8. The equations are set up in such a manner that with increasing n the 
equation with the least value of |^ + 2a?| is in the middle. 
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FIGURE 1: Modulus and sign o/ c§ n (h 2 )/c§(h 2 ) as a function of v for two values of the parameter h 2 , 
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To avoid overflow in applying (3), this relation is rewritten as a recurrence formula for the 
quotients c£ n /c£ w ._ 2 , given by 



cgn+2 ^-(v+2n) 2 / cY ln 

* A 2 V*.-t 



= 0; 



(4) 



compare [7]. The starting values c% N lc£ N _ 2 and c v _ 2S lc v _ 2N+2 were evaluated from continued fraction 
expansions obtained from (4) with c£ (A>m) and c" N+m) set equal to zero. A numerical test deter- 
mines the number of elements, ra, to be used. 

Figure 1 shows that as a function of v, cgjcfi is continuous for noninteger v. It has n zeros; 
they lie near integer values of v for small A 2 (fig. la-c). For increasing values of A 2 the zeros move 
toward smaller v (fig. Id— f). For interger v, however, c%Jcl is discontinuous. The relation 



c?l m (h 2 ) = c? r - m (h 2 ) 



(5) 



in which r is an arbitrary integer, is valid for the limiting value from the right, giving k=d m , say, 
andi>=m, ra = 0, 1, . . ..Similarly, 



cp. m (h 2 )=-cT r . m (h 2 ) 



(6) 
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Figure 2: Modulus and sign of c|„(l)/cg(l) in the region v ~ m. 

a. m = 3, illustrating c v _.Jc^ and c"_Jc^ b. m =4, illustrating r'L.Jc^ and c l '_Jc^ 
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is valid for the limiting value from the left, giving \=b m , say, and v=m, m=l, 2, . . .; compare 
[1], p. 116, and [2]. In particular, for even values of m we have from (6) 

c? m (h?)=^c? m (h*)=0. (6a) 

Fourier coefficients with other indices coinciding according to eqs (5) and (6) with j>— » m are shown 
in figure 2. The absolute values of the Fourier coefficients with negative indices vary rapidly in 
the region m — e < v < m,e <^ 1. 

Tamir and Wang [2] calculated only the zeros of c^Jcq which result from eq (6a), although by 
use of (5) and (6) it is also seen that for v¥"m c\Jc\ can vanish. Thus the graphs in [2] become 
false. This leads to an erroneous calculation of Mathieu functions with indices in the neighborhood 
of integers. 

Summary 

A method for calculating the Fourier coefficients of Mathieu functions in stable regions 
based on the method of bisection and Miller's recurrence algorithm is introduced. Some Fourier 
coefficients are calculated and compared with those given earlier by Tamir and Wang. It is shown 
that the method of Tamir and Wang fails for indices in the neighborhood of integers. 
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